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Abstract. The polarization of radiation by scattering on an atom embedded in 
combined external quadrupole electric and uniform magnetic fields is studied the- 
oretically. Analytic formulae are derived for the scattering phase matrix. Limiting 
cases of scattering under Zeeman effect, and Hanle effect in weak magnetic fields 
are discussed. 
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1. Introduction 

Scattering of polarized radiation by an atom is a topic of considerable interest to as- 
tr ophysics, especial ly with the advent of imaging polarimeter systems like ZIMPOL I & 
II IjGandorfer 20031) which reach accuracies of the order of for measuring the Stokes 
parameters characterising the observed radiation. The concept of scattering phase ma- 
trix connecting the Stokes vector S' of incident radiation to the Stokes vector S of scat- 
tered radiation was introduced quite early in the context of Rayleigh scattering. Landi 

Landi Degl'Innocenti, Bommier and Sahal- 



DegPInnocenti (1983a: 



Brechot (1990; 



1991a; 



1983b; 



1984; 




1991b) and Bommier (1997a; 



1997b) developed a comprehensive 



theoretical framework to describe the generation and transfer of polarized radiation in 
spectral lines, formed in the presence of an ex ternal magnetic fi eld. I n the context of 
radiation transfer work, Stenflo and Stcnholm 1 1976 ) and Rees jl97s| . used co mplet e 
frequency redistribution (CRD) in the resonance scattering, and Dumont et al. (1977), 
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Rees and Saliba p982(l . Nagendra <jl986; 



1988), Faurobcrt (1987), Ivanov et al. (1997) 



and later works employed partial frequency redistribution (PRD) line scattering mech- 
anisms in the absence of magnetic field. The Hanle effect is a depolarizing phenomenon 
which arises due to 'partially overlapping' magnetic substates in the presence of weak 
magnetic fields, where the s plitti ng produced is of the same order as or less than the 
natural widths. Favati et al. (1987) proposed the name 'second Hanle effect' for a similar 



effect in 'electrostatic fields'. Casini and Landi Degl'Innocenti (|1993l) have discussed the 
problem in the presence of electric and magnetic fields for the particular case of hydrogen 
lines. The relative contributions of static external electric fields, motional electric fields 
and m agneti c fields in the case of hydrogen Balmer lines, have been studied by Brillant 
et al. (1993). A historical perspective and exte nsive references to earlier lit eratur e on 
polarized line scattering can be found in Stenflo |l994L Trujillo Bueno et al. I 2002 ) and 



Landi Degl'Innocenti and Landolfi (2004). The purpose of this paper is to derive the 
scattering phase matrix for the case of combined magnetic and electric quadrupolc fields 
with arbitrary strengths. The particular case of transitions between J — and J = 1 
states is considered following Oo et al. (2004; 2005). 



2. Theoretical Formalism 

The Hamiltonian for an atom, when it is exposed to an exter nal magnetic fie ld B together 
with an arbitrary external Coulomb field is of the form IjOo et al. 2004|) 



H = H +gJ-B + J2 {-l) m Vi m Qi- m , (1) 

1=0 m=-l 

where g denotes the gyromagnetic ratio, J the total angular momentum operator for 
the atom with components J x , J y , J z while VJ m denote the 2'-pole components of $ and 
the components Qi- m characterise the electric charge distribution inside the atom. It 
is well-known that the magnetic field B splits an eigenstate of Hq with energy E and 
total angular momentum quantum number J into (2 J + 1) equally spaced levels | J M) 
with energies Em = E + gBM, where B = \B\ and M denotes the magnetic quantum 
number with respect to an axis of quantization chosen along B. If the atomic states 
are eigenstates of parity, the I = 1 term makes no contribution. In an external electric 
quadrupolc field with 1 — 2, Vim may equivalcntly be expressed in terms of the cartesian 
components V ai p with a, (3 = x, y, z of a traceless symmetric second rank tensor, which 
defines its own Principal Axes Frame (PAF), wherein V a .p = V a , a Sa,p-> so that 

2 



Q 2 -rn = A 3 — J ■ J + 7} { J% - J 2 y } 



m=-2 



(2) 



where A — QV ZZ /{4J(J— 1)}, if Q denotes the electric quadrupole moment of the atom 
and -q = (V xx — V yy )/V zz denotes the asymmetry parameter of the field. In such a case, 
the (2J + 1) substates \tpk) where fc = 1, 2, • • • , (2J + 1) with energies Ek are neither 




Fig. 1. The scattering geometry: (Xq,Yq,Zq) refers to the Principal Axes Frame 
(PAF) characterising the electric quadrupole field. The radiation is incident along (#', cf>') 
and scattered along (#, <j)) with respect to the astrophysical reference frame denoted by 
(X, Y, Z). The magnetic field B is oriented along (Qb : <Pb) with reference to PAF and 
{0b,4>b) with reference to the astrophysical reference frame (the azimuthal angles 4>b 
and <f>B are not marked in the figure). 

equally spaced nor are they identifiable as \J M) states. We may, however, represent 
them as 

J 

m= E <i(B,A, v )\JM) , (3) 

M=-J 

in terms of \J M) states defined with respect to PAF, where the expansion coefficients 
a\ { as well as the energies Ek are not only functions of B,A,r] but also of the angles 
(0b,4>b) °f B with respect to PAF (see Fig. P). For a detailed discussion for J = 1, § 
(see Oo et al. 2004; 2005). 

The PAF itself may, in general, be different from the astrophysical frame, in which 
case 

J 

\i>k) = J2 c "l Jm ) ; fc = i,2,--. ,(2j + i) , (4) 

m— — ,7 

in terms of the | J to) states defined with respect to the astrophysical frame and 

J 

<t = Yl a M ( B > A > V) D iu i a Q ,0Q , IQ) , (5) 
M=-J 

if (otQ, Pq,^q) denote the Euler angles of the PAF with respect to the astrophysical 
frame. If the magnetic field alone is present, l^fc) are identical with \ J M) states and 
=Di M (4> B ,9 B ,Q) in Eq. ®. 

We now consider the scattering of polarized radiation by an atom which makes a 
transition from an initial state \tpi) with energy Ei, total angular momentum Jj and 
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parity 7Tj to a final state \ipf) with energy Ef, total angular momentum Jf and parity Ttf 
when polarized radiation with frequency v' is incident on the atom in the astrophysical 
frame in a direction (9',(j)') and gets scattered into a direction (9, <j>) w ith frequrncy v. 



The left and right circular states of polarization as defined by Rose ( 1957J) are denoted by 
fi = ±1. The second order transition matrix element for scattering of polarized radiation 
may then be written, with respect to the polarization states as 

Tw' =E^(l')i4(/i') , (6) 
where the summation is with respect to the intermediate states \il> n) of the atom with 



energy E n , total angular momentum J„ and parity n n . Following l|Oo et al. 2004 ). the 
matrix elements for emission from | J u m u } to | Ji mi) are of the form 

L 

£mim u (v) = ^ C(Ji,L,J u ;m h M L ,m u ) 

L M L =-L 

0^(^,9,0)*, (7) 

where Jl is proportional to the reduced matrix element. The complex conjugate of Eq. lf7| 
defines the matrix elements -Am„m ; (/.t) for absorption of radiation with polarization /x 
incident along [9, (f>) leading to \J U m u ) from | J; mi). Using the notations lu — 2itv 1 lu' — 
2itv'; u> n — E n — Ef and to' n = E n — Ei, the profile function <£>„ = [ui n — u — ir n ) _1 = 
(uj' n — uj' — ir„) _1 where the width associated with \ip n } is denoted by T n and energy 
conservation requires Ei + a/ = Ef + lu. Angular momentum and parity are conserved 
individually during the absorption and the emission. 

3. Scattering Phase Matrix 

If <S' denotes the Stokes vector, which characterizes the state of polarization of the 
incident radiation, the Stokes vector S characterizing the scattered radiation is 

s = n s' , (8) 

where 1Z is a (4 x 4) matrix whose elements are of the form 

K pp ,=J2*nK'K2'> ( 9 ) 

where the phase matrix elements are given by 

?%' = \ E K 7 'W>?W 

fj." ,t* / " 

Tr\£(Li)C n A(p')C A\n"')C n ' £\n")C f ] , (10) 



with p,p' = 0, 1, 2, 3, using the density matrix formalism l|McMaster 1 961) for polariza- 
tion of radiation. The er^ with p — 1,2,3 are Pauli matrices with respect to basis states 
li = ±1 for radiation, while <7q = 1. We use the notation Tr{= J2 m ) to denote the Trace 
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of the (2Jf + 1) x (2Jf + 1) matrix contained within the square bracket. The C\C f ,C n 
and C n are defined in terms of their elements 

pi —J J* . r>f _ J J* 

C n — r n r n * ■ C n ' = r n> r n " ("111 

Each of these matrices C % ,C* ,C n ,C n are clearly hermitian and satisfy the condition 
C 2 = C. Note that the V™, depend not only on n and n' but also on the direction {9b, 4>b) 
and the strength B of the magnetic field B and on (A, rj; ctQ, (3q, 7q) characterising the 
electric quadrupole field (because of Eq. |0J for the with k = i, /, n, n'), apart from the 
angles (8',cf)') of the incident and (9,4>) of the scattered radiation. Explicitly, therefore, 
Vpp,' = Vpp,' {9, <p; 9'(j)'; B; A, 77; a Q , /3q,7q) for any given J t ,J f . In the case of resonance 
scattering, when only a single intermediate state \ip n ) with E n — Ei + lu' = Ef + ui 
contributes to Eq. ©, one can replace the summation over n, n! by n — n 1 corresponding 
to a single excited level, whereas the double summation over n, n' has to be retained in 
Hanle scattering. 



4. Particular Case 



We consider the simple case of scattering with electric dipole transitions between a total 
angular momentum zero lower level and a total angular momentum one upper level, i.e., 
J; = Ji = Jf = Q and J u = J n = 1. Clearly, C % = C s = 1. We may then use Eq. @ to 
simplify the product £* (//') £(/**) and use the complex conjugate of Eq. Q to simplify 
the product A(jJ) A^ (//") in Eq. (JTDJl, so that 



V™: = \gf{9,4>)g n /(8' : 4>y ; 

in terms of 

2 1 

s; n >,0)Hi7i| 2 E 2 C(l,l,X;m',-m",m) 

A— rn/ ,m ;/ — — 1 

A 



n, n' = 1, 2, 3 , 



2 

71 



5\,o + 5\,2 



where 
/o(A,0 

/l(A,£) = -5a, 2 <5?,2 + <Jf,-2 

/2(A, = i <5 A ,2 <5?,2 - <5 5: _2 

/ 3 (A,0 - \/2^, 5 A)1 . 



In the absence of the electric field, n = M, n' = M' and the c^f = D\ 



(12) 



(13) 



(14) 

leading to the we ll-kn own Hanle scattering phase matrix given by Eqs. (9) to (16) of 
Landi and Landi Ijl988|) if <&a/ ®*m> can be assumed to be independent of M, M 1 in the 
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limiting case of weak fields. If the Doppler convolution is effected following exactly the 
procedure outlined by Stenflo (1998), the Hanlc-Zceman scattering matrix represented 
by Eqs. (49) and (50) of Stenflo (1998) is recovered for 9b = <Pb = 0. In the case of 
strong fields i.e, if g B is large compared to the line widths (Zeeman effect), one may set 
M = M' and recover Eq. (52) of Stenflo (1998) for Ob — 4>b — and the results obtained 
much earlier by Obridko |l965aj). 

5. Numerical Results and Discussion 

If we consider the simplest geometry of the combined magnetic and quadrupole electric 
fields with B along the Z-axis of the PAF with the PAF itself coinciding with the astro- 
physical frame i.e., ctQ — [3q — 7q = 6b = <t>B = Ob = 4>b = 0, the upper level with 



J u = 1 is split into three levels n — 1,2,3 l|Oo et al. 20 041 with energies 



E 1 = -2rgB; E 2 , 3 = (r T a)g B , (15) 

where r = A/gB the ratio of the electric quadrupole and magnetic field strengths, s = 
(r 2 r] 2 + l) 1 / 2 and the corresponding eigen states l^i): an d characterized by 

r 1 - 1 r 1 - r 2 ' 3 - 

c — 1, c ±1 — c — u 
4i = ±4l = (s + r7 1 ±l)/2(s 2 + r7 ]S ) 1/2 , (16) 

using Eqs. Q and JSJl, with the electric quadrupole field strength A > 0. 

To understand the combined effect of magnetic and electric quadrupole fields, we 
present in Figs. J3Ji-[2f) the general behavior of the scattered Stokes line profiles S for 
a given unpolarized incident radiation, S = (1 0) T , for particular choices of the 
directions (0' = ?, <\>' — 0) and (0 = f , <j> = tt). We compare these with the pure Zeeman 
scattering case (see Figs. EK"Et) - In the Stokes Q profile, the positive maximum at the 
line center and the negative maximum symmetrically placed at a±\ components, which 
are typical of the well known Zeeman effect. The maximum of the V profile at a±\ 
components have opposite sign, which is also a well known characterisitc of Zeeman 
effect. We assume here that the magnetic field and the quadrupole electric field are 
equally strong (i.e., r = A/B = 1). We also assume B to be four times the natural 
line width and set n = n' = 1,2,3. The solid and dashed curves in Figs. J2i-|2f) are 
computed for the values of 77 = and 1 respectively. In the combined fields case, the line 
component arising due to the \ipi) = |1, 0) state (which represents the central component 
in the corresponding pure Zeeman case) is positioned in the red wing (see Figs. |2i-|2F). 
The unequal strengths of the scattered line profiles arising from \ip2) and ^3) states are 
clearly seen in all the scattered Stokes line profiles (/, Q, V). This is due to the weighted 
superposition of the magnetic substates |1, — 1) and |1, 1). In the scattered Stokes V line 
profile, \^)\) state does not contribute, as in the case of pure Zeeman scattering. Therefore 
the shape of the scattered Stokes V profile is similar to the Zeeman case, except for the 
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Fig. 2. Stokes line profiles in arbitrary units comparing the Zeeman scattering in a pure 
magnetic field (panels a — c) , and scattering under the combined magnetic and electric 
quadrupole fields (panels d — f ) with r = 1. The solid lines correspond to asymmetry 
parameter r\ = 0, and the dashed lines to r\ = 1. In all panels, x is the frequency 
displacement from the line center in natural width units. 



shifting of the position and the change in the relative strength of components when r\ > 0. 
If i] vanishes, the strengths of both components are same. 
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